Draft version August 12, 2011 

Preprint typeset using I^T^X style cmulatcapj v. 11/10/09 



CALIBRATING THE STAR FORMATION RATE ATZ~1 FROM OPTICAL DATA 

Nick Mostek 1 , Alison L. Coil 2 , John Moustakas 2 , Samir Salim 3 , and Benjamin J. Weiner 4 

Draft version August 12, 2011 

ABSTRACT 

We present a star-formation rate calibration based on optical data that is consistent with average 
observed rates in both the red and blue galaxy populations at z ~ 1. The motivation for this study 
is to calculate SFRs for DEEP2 Redshift Survey galaxies in the 0.7 < z < 1.4 redshift range, but 
our results are generally applicable to similar optically-selected galaxy samples without requiring UV 
or IR data. Using SFRs fit from UV/optical SEDs in the AEGIS survey, we explore the behavior of 
restframe B-band magnitude, observed [O II] luminosity, and restframe (U-B) color with SED-fit SFR 
for both red sequence and blue cloud galaxies. The resulting SFR calibration produces fit residual 
errors of 0.3 dex RMS scatter for the full color-independent sample with minimal correlated residual 
error in L[0 II] or stellar mass. We then compare the calibrated z <~ 1 SFRs to two diagnostics 
that use L[0 II] as a tracer in local galaxies and correct for dust extinction at intermediate redshifts 
through either galaxy i?-band luminosity or stellar mass. We find that a L[0 II] - Mb SFR calibration 
commonly used in the literature agrees well with our calculated SFRs after correcting for the average 
£?-band luminosity evolution in L + galaxies. However, we find better agreement with a local L[0 II]- 
based SFR calibration that includes stellar mass to correct for reddening effects, indicating that stellar 
mass is a better tracer of dust extinction for all galaxy types and less affected by systematic evolution 
than galaxy luminosity from z = 1 to the current epoch. 

Subject headings: galaxies: evolution — galaxies: active — galaxies: high-redshift 



1. INTRODUCTION 

The global star formation rate (SFR) of large galaxy 
samples at intermediate redshi ft has been an area of ac- 
tive study for the past decade (pteidel et al.|fl 999; Hop- 
kins fc Beacoinl[2006l |Zhu, Moustakas fc Blanton]|2009| 
Gilbank et al.||2010* ). Each SFR study has chosen a par- 
ticular star-formation rate indicator (or combination of 
indicators) from a variety of possible measures. The in- 
dicators associated with star formation cover a broad 
range of wavelengths of the electromagnetic spectrum 
from x-ray to radio. The most studied measures include 
UV luminosity, which indicates the amount of massive 
star formati on within the compos ite stellar population 



of a galaxy ( |Heckman et al.||1998 [), and the luminosity 
of nebular emission lines such as Ha and [O II], which 
measu re the amount o f ionizing radiation f r om massive 
stars dKennicuttl [19981 iKewley et al.||2004| |Moustakas, 
Kennicut fe Tremonti||z006[ hearafter MU6). Addition- 
ally, deep observations from Spitzer have tied the mid-IR 
continuum emission at 2 4/^m to compact star formation 
regions in l ocal galaxies ( Calzetti et al.|[2~007 Kennicutt 
et al.||2009[ ). 

AH SFR diagnostics essentially measure the energy 
output of young stars either through direct observation 
of the UV flux or indirect observation of the ionizing lu- 
minosity through nebular emission lines or reprocessed 
emission from dust. Regardless of the diagnostic used, 
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the SFR indicators have systematic effects which must 
be removed in order to give an accurate measure of the 
SFR. The primary physical effect that dominates system- 
atic uncertainty in restframe UV/optical observations is 
the amount of dust extinction in the host galaxy. Galaxy 
dust extinction attenuates the UV continuum or emission 
line luminosity and biases the indicators to lower SFRs 
I KennicuttJ|1998[ jKewley, Geller fc Jansen||2004] |Kew 



ley et al.||2004[ ). Further, the amount of dust extinction 



within the star-forming gala xy population also co rrelates 
with galaxy luminosity (Wang & Heckma n|1996|), stellar 
mass jBrinchmann et ai. 2004; N OLvsk e^FaT^UUTb] Elbaz 
et al. 2007; Gilbank et al. 20I0|), and star- formation his- 
tory ( JNoeske et al.||2007a[ ), presumably due to increased 
metallicity and dust production from massive stars and 
mass-loss return to their ISMs. Producing accurate dust- 
corrected SFRs from UV continuum or nebular emission 
lines requires at minimum a measure of the dust extinc- 
tion through diagnostics such as the UV/optical slope 



(Calzetti et al. 1994 Kong 2004) or the Balmer decre 



ment (Ha/H/j) assuming case B recomb ination for typi 



cal galaxy environments ( Calzetti]|2001 1. 

Additional systematic uncertainties can also limit the 
SFR accuracy after dust extinction has been removed. 
For example, emission line luminosities from [O II] and 
[O III] depend on other nebular characteristics such as 
the chemical abundance and excitation states of the gas 
(Kewley & Dopita 2002). M06 estimates that even af- 
ter an ideal Balmer-decrement dust correction is applied, 
[O II]-based SFRs have a lower error limit of 32% for 
galaxies with typical metallicities and can have consider- 
ably larger error for more extreme abundances. Further, 
the gas can be excited from sources that are not associ- 
ated with actual star formation, such as the hard-ionizing 
background radiation created by AGN, which can be a 



2 



dominant effect in quiescent galaxies ( Yan et aL||2006 ). 

Beyond the UV and optical re gime, IK measurements 
can be used to trace galaxy SFR ( Kennicutt|19 98 ; Daddi 



et al.|2007 ); UV photons from massive stars are absorbed 



by dust within the host galaxy ISM and re-emmited 
at IR wavelengths. IR-based SFRs are most accurate 
for heavily attenuated environments such as young star- 
bursts where dust-corrections and measurements of the 
UV luminosity are difficult. However, IR observations 
tend to underestimate the dust luminosity in normal spi- 
rals where much of the UV radiation is unabsorbed. IR 
luminosities can also have systematic variations due to 
differing dust geometries in the host galaxy or contami- 
nation from old s tellar populations , particularly in red- 
sequence galaxies ( Sal^er^]2009 hereafter S09) . Nev- 
ertheless, recent work has calibrated the SFR to I R lumi- 



no sity in various bands for these dusty g alaxies ( Zhu et 
aT|2008| |Calzetti|2009| |Rieke et al.||2009| ) and calibrated 



the 24/xm luminosity to Ha luminosity to provide a mea- 
sure of both the SFR and account for typical dust atten- 
uation encountered acr oss a broad range of blue galaxies 
(iKennicutt et al. 112009). 



One way to mitigate these various systematic effects 
is to obtain integrated spectral and photometric mea- 
surements over a wide wavelength range, providing con- 
straints on the dust extinction and independent cross- 
checks in SFR. In S09, galaxy template SEDs were fit 
to measurements spanning FUV to iVband wavelengths 
available from the All- Waveleng th Extended Groth Strip 
International Survey (AEGIS) ( |Davis et aL||2007| ). The 
S09 SFRs were based primarily on measurements of the 
UV luminosity and are combined with a Baysian analysis 
of stellar population synthesis models to fit the template 
SEDs. By simultaneously fitting across the panchro- 
matic data and correcting for various astrophysical ef- 
fects such as dust extinction, UV upturn in red galax- 
ies, and SFR timescales, the dust-corrected SFRs are 
robust against large systematic errors that would typi- 
cally limit SFR measurements at intermediate redshifts 
across galaxy types. We consider the SED SFRs fit to 
UV/optical measurements to be the least biased mea- 
sure of SFR in red galaxies, particularly when the UV 



Han et al.|[2007 ). These 
:erence for the SFR cali- 



upturn is taken into account 
SED-based SFRs provide a re 
bration presented here and will be discussed further in 
Section GO 

The goal of this work is to take advantage of these 
robust SED-fit SFRs and provide a SFR calibration 
for optically-selected galaxy samples where wide multi- 
wavelength coverage or emission line data may not be 
available. Specifically, we wish to develop a suitable cal- 
ibration using SFR data in the AEGIS field and extend 
that calibration to z ~ 1 galaxies in the DEEP2 Redshift 
Survey (Davis, 2003) outside the AEGIS field. Because 
the SFR calibration is based on a large galaxy sample, 
the final means-tested SFR calibration will not be highly 
accurate for individual galaxies. However, it will repro- 
duce the global SFR trends seen in the z ~ 1 AEGIS 
sample and allow us to confidently study the galaxy envi- 
ronment at these redshifts with respect to SFR. In <p]we 
describe the volume-limited galaxy sample selected from 
the DEEP2 survey that is matched to the S09 AEGIS 
SFRs. In Sj3] we investigate which combination of DEEP2 



band luminosity, or restframe color, that are most corre- 
lated with SFR and will therefore provide the best lever- 
age in calibration. Our SFR model for these tests uses a 
weighted linear combination of observed parameters (not 
corrected for dust extinction), and we investigate which 
parameter combination delivers the most accurate SFR 
calibration for gal axie s wit h bo th red and blue optical 
restframe colors (p^ijand p3| . In jg) we describe two 
empirical SFR estimators which use local [O II] luminos- 
ity measurements and correct the SFRs for systematic 
effects via Mg or stellar mass. We then compare these 
diagnostics to our best-fit SFR calibration at interme- 
diate redshifts, and we present our conclusions in Sj5| 
Thr oughout this wo rk, we adjust all masses and SFRs 
to a Salpeter ( 1955 ) IMF and assume a ACDM cosmol- 
ogy (SJm = 0.3, J 2a = 0.7, h = 0.7). All magnitudes used 
in this study are in the AB system. 

2. DATASETS 

The main sa mple for this work is drawn from the 



2003) where the [O II] 
been measured with the 



DEEP 2 survey ( |Davis et a! 
doublet emission lines have 
DEIMOS spectrograph on Keck ( |Faber et al.|2003| . The 
DEEP2 survey contains four fields that are separated 
across the Northern sky for year-round observation, and 
one of the fields (Field 1) overlaps the AEGIS survey 
in the Extended Groth Strip (EGS). Within the EGS, 
DEEP2 targets all galaxies to Rab < 24.1 and is limited 
to z < 1.4, beyond which the [O II] emission line falls out- 
side the observed wavelength range (6500-9200A). In the 
non-EGS fields, targets selected by DEEP2 are limited 
by design at a lower redshift of z > 0.7 due to an opti- 
cal BRI color cut. The total DEEP2 sample with good 
redshifts (> 95% confidence, z quality flag > 3) is 31,656 
galaxies. When available, the emission line equivalent 
widths in DEEP2 are measured using a nonlinear least- 
squares fit of a Gaussian to the emission line profiles. A 
line flux is then computed by measuring the continuum in 
a 20-60A window around the line and calibrating the flux 
to iVcorrected RI photometry. This method takes into 
account both the throughput and slit losses for DEEP2 
galaxies of which nearly all have an effective radius less 
than the 1" DEIMOS slit width (We iner et. al, 2012, in 
preparation; description of metho d in Wcincr et al.|2007 
Zhu, Moustakas fe Blantonfl2009[ ). 

To build a sample of z ~ 1 SFRs, we match the DEEP2 
galaxy sample to available panchromatic data from the 
AEGIS survey. The matched AEGIS sample contains a 
total of 5345 objects and is a subset of the total avail- 
able DEEP2 spectra in the EGS field. From this sample, 
we obtained SFR estimates from template SEDs which 
are fit to UV, optical, and IR -FT-band photometry in 
S09. The composite templates are c onstructed from stel- 
lar po pulation synthesis models of |Bruzual fe Chariot] 
(20031 augme nted with a two-comp onent dust attenua- 



tion model of Chariot fc Fall] ( 2000 1 to account for dust 
associated with various star-formation histories and the 
intergalactic medium. The templates used in S09 also 
take into account extreme blue horizontal branch stars 
which are thought to be responsible for the "UV up- 
turn" seen in s ome early-type ga laxies with old stellar 
populations (see Salim et al. 12007 for details of the SED 



models). For 45% of the S09 sample, the SED-fit SFRs 



measurements, observed [O II] luminosity, restframe B- have been compared with 24/im measurements for both 
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red sequence and blue cloud galaxies. The total IR lu- 
minosities are inferred from the 24/im flux densitie s and 
galaxy redshift s by fitting to IR SED templates (Dale 
fc Helou|[2002| ). S09 found that SED SFRs of theTTue 
galaxy population (NUV-R< 3) were consistent with the 
measured IR luminosities to 0.3 dex RMS. The SFRs 
for red galaxies (NUV-R> 5) were less well matched in 
IR luminosity due to the increasing contribution of light 
from old stellar populations in the IR measurements. Be- 
cause S09 accounts for the latter systematic effect though 
template fitting, the S09 SFRs based on UV SEDs are a 
better estimate of the true SFR than the 24/xm measure- 
ments in red galaxies. 

The stellar masses for our gal axy sample are com- 
puted from the Bell et aT| ( 2003 ) color- M/L relations, 
modifi ed to DEEF2 redshifts through the |Weiner et al.| 
(2009) calibration. The color-M/L calibration is a lin- 
ear set of equations that require several measured quan- 
tities, include Mb, restframe (U-B) and (B-V) colors, 
and galaxy redshift. Independently, we confirm that 
our calculated stellar masses from these relations agree 
with stellar masses derived from if -band measurements 
( |Bundy et aTp 006) to 0.3 dex RMS with no systematic 
offset. As an additional consistency check, we find that 
the color-M/L stellar masses agree with those generated 
by the more sophisticated S09 SED-fitting techniques to 
within 0.3 dex. We use stellar masses from the M/L-color 
calibration rather than stellar masses obtained from the 
SED fitting because we want the input parameters in our 
SFR calibration to be easily reproducible and indepen- 
dent from the SED fitting that yields our fiducial SFRs. 

2.1. Sample selection 

To produce a statistical sample volume-limited in 
for both red and blue galaxies, we restrict the matched 
DEEP2 / AEGIS galaxy sample to M B < -20 for 
0.74 < z < 1.0, resulting in 1039 total galaxies. The se- 
lected redshift range and galaxy luminosity limit is sim- 
ilar to color-independent, volume-limited D EEP2 sam- 
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blue galaxies in our sample have no measured amount of 
L[0 II]. The L[0 II] measurement failures in blue galax- 
ies are likely due to interactions of the emission line with 
bright sky emission lines, which happens at an nearly 
equal rate for both blue and red galaxies over this redshift 
range. The incompleteness due to sky line interaction is 
generally independent of the galaxy restframe S-band 
luminosity, color, and redshift. Removing the galaxies 
without measured [O II] reduces the sample size to 885 
DEEP2 galaxies matched with AEGIS and restricts our 
SFR calibration analysis to galaxies with a SED-fit SFR 
limit of log(^)> —1.0 and [O II] line luminosity limit of 
log(L[0 II] )> 39.7 at z ~ lQ 

While blue galaxies in our sample generally have well 
detected L[0 II] measurements (> 3a), a significant frac- 
tion of the red galaxies with L[0 II] measurements have 
large error due to inherently lower line luminosities. Re- 
stricting the L[0 II] measurements to better than 3cr in 
DEEP2 reduces the sample to 80% of the volume- limited 
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Fig. 1. — Comparison of the S09 SFRs for blue galaxies from 
0.74 < x < 1.4 with observed L[0 II] in DEEP2. The plotted points 
are for both uncorrected (black triangles) and dust corrected (red 
diamonds) SFR values. The green stars denote the rejected S09 
SFRs for heavily dust-extincted blue galaxies that require a factor 
> 20 in dust correction. Note that the uncorrected SFRs are tightly 
correlated with observed L[0 II] , indicating that the dust correction 
introduces systematic uncertainty into the corrected SFR values. 



blue galaxies and 27% of red galaxies. For comparison, 



J The SFR ip is in units of Mq yr 
throughout this study. 



, and L[0 II] is in ergs s 



Yan et al. ( 2006]) found that 35% of all volume- limited 
red galaxies in SDSS had [O II] detected at a 3er level in 
EW. Because we are interested in testing SFR calibra- 
tions for a uniform sample with both L[0 II] and Mb, 
we wish not to limit the completeness in Mg, which is 
well measured for all galaxies, due to the higher L[0 II] 
errors. Therefore, we opt to keep all galaxies with mea- 
sured L[0 II] in our sample and weight the SFR fits by 
the measurement errors. 

We place one final restriction on the sample to remove 
galaxies that have a large dust-corrected SFRs in the 
S09 data. Figure [l] shows the dust-corrected and un- 
corrected S09 SFRs in the matched DEEP2 / AEGIS 
sample as a function of observed [O II] luminosity from 
DEEP2. We find that the uncorrected UV/optical SED 
SFRs and observed L[0 II] are tightly correlated due to 
similar dust extinction properties in the UV. After the 
dust correction is applied, a small population (13%) of 
galaxies have SFRs that arc > 3er outliers from the mean 
dust-corrected distribution. These "dusty" star-forming 
galaxies have SFR corrections of more than a factor of 
20 and are noted in S09 as primarily intermediate color 
(green valley) galaxies. Since we are attempting to cal- 
ibrate SFRs at z ~ 1 with limited optical observations, 
we assume that the dust extinction for these obscured 
galaxies cannot be accurately measured on an individual 
galaxy basis (usually done through IR measurements) 
and remove them from the following analysis. While 
removing the outliers allows our calibrations to more 
closely track the "typical" galaxy in our sample, this un- 
derlying assumption restricts the SFR calibrations in this 
study to galaxy samples with moderate dust extinction 
and we acknowledge that our results will not be accu- 
rate for galaxies that require extraordinarily large dust 
corrections. The sample used in our SFR calibration is 
reduced to a total of 771 galaxies after we remove these 
galaxies. 
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Fig. 2. — Histograms of the selected and matched AEGIS / 
DEEP2 sample from 0.74 < z < 1.0 (red, dashed) and the en- 
tire selected DEEP2 galaxy sample over the same redshift range 
in the EGS (black, solid). The DEEP2 EGS sample histograms 
have been normalized to the total number in the AEGIS matched 
sample to show the relative shape of the distributions. The dotted 
vertical lines mark the median values for each sample. 

To confirm that our selected and matched AEGIS / 
DEEP2 sample is a representative subset of optically 
selected z ~ 1 galaxies, we compare the measured pa- 
rameter distributions for the matched AEGIS / DEEP2 
sample to the full EGS DEEP2 galaxy sample using the 
same selection cuts. This test must be restricted to the 
EGS field because the original DEEP2 target selection 
within the EGS was color-independent, whereas the three 
other DEEP2 fields were color selected. Figure [2] shows 
the selected sample histograms for [O II] line luminos- 
ity, restframe B-band luminosity, stellar mass, and red- 
shift. Because the matched and selected galaxy sample 
is roughly 10% of the entire DEEP2 EGS sample, the 
DEEP2 histograms have been normalized to the total 
matched AEGIS sample for easy comparison of the dis- 
tributions. We find that the AEGIS / DEEP2 matched 
sample is an accurate subset of volume-limited, optically- 
selected galaxies in this redshift range and produces 
equivalent median values as the full EGS DEEP2 sample 
with log(L[0 II])=41.5, M B = -20.8, log(M*/M )=lO.3 
and redshift z = 0.88. 

3. SFR CORRELATIONS 

3.1. Trends in blue and red galaxies 

In this study, we explore the systematic correlation be- 
tween the dust-corrected SFR and three measured physi- 
cal parameters: 1) observed [O II] line luminosity, 2) rest- 
frame B magnitude, and 3) restframe (U-B) color. These 
parameters are based only on optical measurements (i.e. 
no UV or IR data iss required) and therefore may be 
readily produced over wide fields where spectroscopic or 
photometric redshifts exist. The correlated trends of ob 
served L[0 II] and Mb with SFR are shown in Figure [3 
where we have separated the least-squares fits of bot 



the measurement error in both the independent and de- 
pendent variables plus an additional " intri nsic scatter" 
term added in quadrature (see Section 3.2 for more de- 
tail). In the case of L[0 II], we plot the measurements 
with < 2tr significance using an open circle to indicate 
that they are less important in the fit. 

Several trends are apparent in the figure: i) (U-B) 
color clearly separates the red galaxies with lower SFRs 
from the blue galaxies with higher SFRs, ii) the red 
galaxies clearly have a greater scatter about the mean 
regression line in both observed L[0 II] and Mb relative 
to the blue galaxies, and Hi) the red and blue galax- 
ies have similar SFR slopes in Mb but have different 
slopes in L[0 II]. The similar SFR slope between red and 
blue galaxies and restframe Mb might indicate that the 
galaxy S-band luminosity is a more accurate indicator 
of the corrected SFR regardless of galaxy color (albeit 
with a small offset between SFR zeropoints), but the 
SFR variation among the red galaxy population is too 
large for a definitive statement. Figure [3] also shows that 
there is a large variation in SFR for blue galaxies of a 
fixed L[0 II], demonstrating that observed [O II] line lu- 
minosities with small measurement errors may still have 
large systematic variations in corrected SFR due to the 
extinction effects. 

The difference in the SFR scatter of red and blue galax- 
ies relative to their mean trends in Mb and L[0 II] occurs 
for several reasons. Red galaxies are typically quiescent 
in star formation and have inherently lower UV flux than 
that of the blue galaxies, leading to higher measurement 
uncertainty in red galaxy SED SFRs. Similarly, the mea- 
surements in L[0 II] have higher measurement error in 
red galaxies simply because the emission line fluxes are 
fainter than blue galaxies. For red galaxies, we also ex- 
pect the measurements are bias ed by [O II ] emission not 



associated with star formation. Yan et al. (2006) showed 
that L[0 II] is dominated by AGlN emission for a signifi- 
cant fraction of the local red galaxy population and that 
this emission is not representative of true star formation. 
We will study the extent to which the observed L[0 II] 
can be used as a direct proxy for quiescent galaxy SFR 
in the following section. 

3.2. Fitting Methodology 

To obtain an empirical SFR calibration for our 
matched DEEP2 / AEGIS sample, we perform a 
weighted \ 2 n t using observed L[0 II] and M^, where 
the galaxies have been separated into red and blue sub- 
samples by restframe (U-B) color. The galaxy SFR fits 
follow a simple linear combination model according to 



n 

Ipmodel = Co + CjPj, 
i=l 



(1) 



where n is the number of observed parameters for the 
linear model and Ci is the scalar fit coefficient for the 
measured parameter p. The linear SFR model is then 
fit to the corrected S09 SFRs, i^swi by minimizing 



fc=i 



(2) 
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Fig. 3. — Correlated trends in SFR with respect to observed L[0 II] and Mj for the restframe red (red circles) and blue (black triangle) 
galaxies in our matched sample. The SFRs used in this study are computed from UV-optical data fit to SED templates in S09 and the 
linear trends are weighted by the measurement errors. The points with error bars in the upper left of each plot show the median error for 
the color-separated galaxy measurements. In the case of L[0 II], open symbols mark the galaxies where the [O II] luminosity has less than 
a 2<r measurement in DEEP2. 



where k is the index of m sample galaxies. The % 2 is 
weighted by the quadrature sum of the observed param- 
eter variance, a 2 , and the intrinsic scatter in the distri- 
bution, of. Typically, the error in ipso9 dominates the 
measured variance when compared to the observed in- 
dependent variables of Mb and L[0 II]. However, in the 
case of the low [O II] luminosity in red galaxies, it is im- 
portant to include the luminosity error to properly weight 
the fit. The intrinsic scatter is computed through a boot- 
strap method by varying the Oi term until the reduced 
X 2 /NDF = 1, ensuring that the computed fit coefficient 
errors are properly scaled and that the fit is unbiased 
given differing errors in the dependent and inde pendent 
variab l es. W e refer the reader to the appendix of |Weiner| 
et al. (2006) for details of weighting fits with intrinsic 



scatter. 

Once the x 2 fit is minimized, we subtract the best-fit 
SFR from the matched S09 SFR and measure the RMS 
of the residual errors. We also judge the performance 
of the calibration by measuring the correlation between 
residual SFR errors and the measured parameters of Mb, 
L[0 II], and stellar mass. A large residual correlation is 
an indication that the calibration process has not fully 
characterized the SFR distribution and that a better cal- 
ibration may be sought. Because stellar mass and rest- 
frame colors are correlated by color-M/L relation, we ex- 
pect that any residual systematic correlated error in color 
from our SFR calibration will also produce a residual er- 
ror correlated with the calculated stellar masses. 

3.3. Color-separated Empirical SFR Calibration 
As discussed in Section 3.1 and shown in Figure [31 we 



initially separate the matched sample by restframeblue 
and red color and fit the linear model of Equation [T] with 
L[0 II] and Mb as the independent variables. We report 
the best-fit coefficients and parameter errors of this cal- 
ibration in Table [l] In Figure [4j we show that by fitting 
only L[0 II] as an independent parameter in each galaxy 



type, we can achieve a SFR fit that produces residual 
errors with an RMS of 0.27 dex scatter for blue galaxies, 
0.46 dex for red galaxies, and 0.31 dex for all galaxies in 
the sample. This result suggests that if restframe colors 
and [O II] measurements are available, L[0 II] can be a 
decent tracer of the average SFR in large galaxy samples 
even though AGN activity may be contributing to the 
[O II] flux in individual red galaxies. However, there re- 
mains a residual systematic correlation along Mb after 
calibration, as indicated by the Pearson correlation coef- 
ficient of r 2 = 0.18. The residual trend indicates that the 
full behavior of SFR has not been calibrated with L[0 II] 
alone and that some improvement can still be made to 
reduce the scatter in the calibrated SFRs. This resid- 
ual correlation is likely due to the degenerate effects of 
dust reddening on the observed [O II] luminosity, B-band 
galaxy luminosity, and restframe color. 

Similarly, we obtain a SFR calibration by fitting the 
linear trends in Mb separately for red and blue galaxies 
(top plot of Fig[5| and the fit coefficients are tabulated in 
Tablcm The blue galaxy SFRs are highly correlated with 
B-band luminosity (Fig. |3| , which provides greater lever- 
age in the calibration over L[0 II]. The fit SFR residual 
errors from the Mb fit show a reduced scatter for blue 
galaxies (0.21 dex), an increased scatter for red galaxies 
(0.54 dex), and no change in the scatter for the combined 
sample (0.30 dex). There is less correlation in Mb (by 
construction) and stellar mass, but we observe that the 
SFR residual errors are now correlated with L[0 II], par- 
ticularly for the red galaxies. Again, the residual corre- 
lations in L[0 II] after calibrating with Mb likely results 
from the differing amount of dust extinction inherent in 
these two observed, uncorrected measures. 

We can also combine both Mb and L[0 II] measure- 
ments to reduce the residual SFR correlations below 
what can be obtained by fitting each independent pa- 
rameter alone. The combined calibration fit achieves an 
RMS scatter of 0.21, 0.45, and 0.27 dex for the blue, red, 
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Fig. 4. — SFR fit residuals using L[0 II] as the independent parameter in the SFR calibration. Fits for red (circles) and blue (triangles) 
galaxies are derived separately. The points with error bars in the upper left of each plot shows the median errors for each galaxy color, 
and the Pearson correlation coefficient values for the combined sample are shown in the upper right. In the case of stellar mass, the error 
on log(M„/MQ) is calculated from the RMS of the DEEP2 masses to the fit S09 masses and therefore accounts for the systematic error in 
estimating stellar mass from the color-M/L relation. 



TABLE 1 

Summary of best fit coefficients for SFR calibration separated by galaxy color 



Fit Parameters 




Blue Galaxies 






Red Galaxies 


Pi 


P 2 


Co, b 




C2,b 


Co,r 




L[0 II] Q 




0.722 ± 0.028 


0.470 ± 0.038 




0.351 ± 0.040 


0.763 ± 0.077 


M B b 




1.099 ± 0.008 


-0.356 ± 0.013 




0.065 ± 0.064 


-0.304 ± 0.083 


L[0 II] 




1.003 ± 0.024 


0.143 ± 0.031 


-0.346 ± 0.014 


0.302 ± 0.055 


0.738 ± 0.081 -0.088 ± 0.069 



a The [O II] luminosity has a zeropoint of L[0 Il]=10 41 ergs s 1 cm 
b The B-band luminosity has a zeropoint of Ms =-21. 



TABLE 2 

Intrinsic scatter and residual errors for SFR 
calibrations separated by galaxy color 11 



Parameters 


Intrinsic Scatter 




Residual RMS 


Pi 


Pi 


blue 


red 


all 


blue 


red all 


L[0 II] 




0.22 


0.29 


0.00 


0.27 


0.46 0.31 


M s 




0.15 


0.49 


-0.01 


0.21 


0.54 0.30 


L[0 II] 


Ms 


0.13 


0.29 


-0.01 


0.21 


0.45 0.27 



a SFR errors are in units of dex. 

and combined galaxy sample, respectively (c.f. bottom 
of Fig. [5]). The inclusion of L[0 II] with Mb actually pro- 
vides constraint on the red galaxies, reducing the scatter 
from 0.54 to 0.45 dex. This result may seem counterin- 
tuitive since the measurement error of L[0 II] is larger 
for red galaxies than blue galaxies and AGN are known 
to contribute significantly the [O II] luminosity. How- 
ever, it is worth noting that galaxies with red restframc 
colors will span larger range of dust extinction values 
and SFR, from elliptical galaxies with low dust content 
and low dust-corrected SFR values to late-type galax- 
ies with heavy dust extinction and high dust-corrected 
SFR. Even with significant measurement error, we find 
that the additional information from L[0 II] to Mb helps 
provide a valuable constraint on the dust-corrected SFR 
in red galaxies. The best fit coefficients for all three color- 



separated SFR calibrations presented in this section are 
summarized in Tableland the residual RMS scatter and 
mean residual offsets are reported in Table [2j 

3.4. Constructing a single SFR calibration 

While separating the red and blue galaxy populations 
enables a SFR calibration across galaxy types with min- 
imal residual correlation and scatter, it inherently re- 
quires all three measurement parameters considered in 
this study (L[0 II] , Mb, and restframe color) to utilize 
the relation. Ideally, one would like a single, simple re- 
lation to calibrate all galaxies to the SED SFRs while 
minimizing the number of required measurement param- 
eters. As shown in Section pO| the red and blue galaxy 
SFRs clearly have differing trends in the observed pa- 
rameters and the intrinsic color dependence cannot be 
ignored. For example, a SFR calibration using L[0 II] for 
all galaxies without color information produces a 35% in- 
crease in the SFR fit residual scatter with large residual 
correlations in L[0 II], Mb, and M*. As an alternative to 
separating galaxies by restframe color, we can add rest- 
frame color information directly into the SFR calibration 
and fit the entire galaxy sample simultaneously. Includ- 
ing the restframe color as an independent parameter in 
the fit produces a continuous SFR calibration between 
galaxy types and helps break the degeneracies between 
the galaxy color, SFR, and dust reddening. 

Figure [3 demonstrated that the dust-corrected SED 
SFRs are "highly correlated with galaxy B-band hum- 
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Fig. 5. — SFR fit residuals using only (top) and both L[0 II] and Mj (bottom) as the independent parameters in the SFR calibration. 
Fits for red and blue galaxies are derived separately. 

nosity, which provides excellent calibration leverage for 
blue galaxy SFRs and marginal leverage for red galaxy 
SFRs. We also observed that the correlation between 
SFR and Mb between red and blue galaxies may have 
similar slopes with only a difference in the zeropoint of 
the relation. Motivated by these observed trends, we 
use the form of Eq. [I] and perform a fit with both rest- 
frame Mb and (U-B) color for all galaxies; the resulting 
residual errors are shown in Figure M (top plot). We 
find that the SFR behavior of the fulTgalaxy sample is 
not well described by a simple linear fit in Mb and (U- 
B) alone. The SFR calibration produces residual errors 
with a a u = 0.33 dex scatter for the full sample and large 
residual error correlated in Ms and stellar mass. The 
residual correlations are apparently non-linear, particu- 
larly in the case of stellar mass, and therefore are not 
fully characterized by the Pearson correlation coefficient 
alone. 

Because the fit residuals in the and (U-B) color cal- 
ibration are seemingly non-linear, we modify our linear 
SFR calibration model (Eq. [IJ to include a second-order 
term in (U-B) (bottom plot of Figure [61). The inclusion 
of this non-linear color term produces an better fit to the 
dust-corrected SED SFRs for the full galaxy sample. The 
final RMS scatter in the SFR fit residuals are (Jbiue = 0.19 
dex, u re d = 0.47 dex, and a a ii = 0.26 dex. The fit resid- 
uals from this calibration are statistically better than all 
other methods considered in this study, and the calibra- 
tion does not require L[0 II] to be included in the fit 



as was necessary in |3.3[ Further, the SFR calibration 
leaves no significant correlations in Mb, M*, and L[0 II] 
(r 2 < 0.032). 

The reason for the poor fit using only linear terms in 
Mb and (U-B) alone can be seen by studying the behav- 
ior of the galaxy SFRs in the Mb - (U-B) color plane. 
Figure [7] shows the S09 SFRs matched to sources within 
our DEEP2 sample from 0.74 < z < 1.4; the overplot- 
ted contours correspond to lines of average SFR with 
log(</0=[0.3, 0.6, 0.9, 1.2] M yr" 1 . The SED-fit SFR 
behavior are decidedly non-linear in Mb and (U-B) color. 
For the blue galaxies, the SFRs primarily increase with 
Mb with only a weak dependence on (U-B) color. This 
SFR behavior in the blu e cloud is in agreeme nt with mea- 
sured SFRs at z = 1.4 (iWeiner et al.||2009|) and is con- 



sistent with the assumption that the U V s 



ope and (U- 



B) colors are correlated, therefore associating Mb with 
UV flux. However, as star formation trends from active 
to quiescent around (U-B)=l, the SFR contours rapidly 
transition to much lower SFRs and have a weaker depen- 
dence on Mb. 

In addition to the differing SFR trends in Mb between 
red and blue galaxies, the differential extinction between 
independent fit parameters, such as the uncorrected B- 
band luminosity and observed [O II] luminosity, can cou- 
ple the fit parameters such that independent linear co- 
efficients are not an accurate description of the global 
SFR trends. Using an additional second-order (U-B) 
color term in the SFR calibration allows an extra degree 
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Fig. 6. — SFR fit residuals using Mb and (U-B) as the dependent parameter in the SFR calibration for all galaxies. The top figure shows 
the fit residuals from a linear fit in both independent variables, while the bottom figure shows the result of a second-order polynomial in 
(U-B). 

ure[8]). The fit residuals in this "multicolor" calibration 
are very similar to the second-order (U-B) calibration, 
with a a u = 0.27 dex scatter and small residual correla- 
tions in all three measurement parameters (r 2 < 0.045). 
Again, we see that an additional color parameter allows 
the model more freedom to decouple the intrinsic galaxy 
color and the effects of dust in the SFR. The combined 
sample SFR calibration coefficients and fit residual errors 
for each of these fit models is tabulated in Table [3] 

We note that it is not entirely surprising that a useful 
empirical SFR calibration can be achieved using observed 
B-band luminosity and multiple restframe colors; these 
measures are similar to the quantities that construct the 
SED-fit SFRs. Either the multicolor or second-order (U- 
B) calibration can be used to produce statistically similar 
SFR calibrations within our volume-limited (0.74 < z < 
1) DEEP2 / AEGIS sample. However, we will see in 
the following section the multicolor calibration is better 
behaved in the Mb - (U-B) plane when extended to the 
full 0.74 < z < 1.4 DEEP2 sample. 

3.5. Trends in Mb - (U-B) plane 

An additional comparison diagnostic of our SFR cali- 
bration models can be observed in the Mb - (U-B) plane. 
Figure [9] shows the SFRs calibrated from Mg and L[0 II] 
using the bimodal restframe color-selected calibration de- 
tailed in Section [373) The galaxy data points are taken 
for the full DEEP2 sample from 0.74 < z < 1.4 and 
are color-coded with Alog(VO=0.3 dex SFR bins. The 




Fig. 7. — The M B - (U-B) plane with SFRs measured in the 
AEGIS survey (S09) and matched to the DEEP2 catalog from 
0.74 < z < 1.4. The contour lines correspond to average SFRs 
of log(VO=[0.3, 0.6, 0.9, 1.2] M yr" 1 . 

of freedom within the fit to simultaneously correct for 
both intrinsic galaxy color between red and blue galaxies 
and dust reddening effects. By separating the degener- 
ate color effects with an additional color fit parameter, 
we can effectively reduce the systematic coupling in the 
SFR fit residuals. 

Knowing that additional color information can produce 
an accurate SFR calibration, we also generate a fit to the 
SED SFRs using M B , (U-B), and (B-V) colors (see Fig- 
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Fig. 8. — SFR fit residuals using MB, (U-B), and (B-V) as the independent parameters in the SFR calibration for all galaxies. The 
residual errors are tabulated in Table \3\ 
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Fig. 9. — The diagnostic M s - (U-B) plane with SFRs calibrated 
using the - L[0 II] calibration split between red and blue galax- 
ies (dotted line). The solid lines correspond to constant contours 
of SFR generated from the matched S09 AEGIS galaxy sample 
with log(i/>)=[0.3, 0.6, 0.9, 1.2] M yr" 1 and are plotted here for 
easy reference to the original calibration data. See Section |3.3| for 
details of the color-split SFR calibration. 

overplotted solid lines correspond to contours of constant 
SFR with log(V>) = [0.3, 0.6, 0.9 1.2] M yr^ 1 from the 
S09 AEGIS sample (c.f. Figure j7l). In general, the color- 
separated SFR calibration tracts the S09 contours well, 
particularly in the blue galaxies where the SFR depen- 
dence is strongly correlated with Mb and independent of 
rcstframe color. However, the drastic split between the 
galaxy colors near the green valley is clearly not represen- 
tative of the true SFR transition at (U-B)~l. Addition- 
ally, bright galaxies just below the (U-B)=l line have 
overpredicted SFRs relative to the S09 contours while 
bright red galaxies just redward of the line have under- 
predicted SFRs. Figure [7] shows that there is some "mix- 
ing" of SFRs for the brightest galaxies with optically red 
rcstframe colors, indicating that galaxies with (U-B)> 1 
at these redshifts can either be old early-type galaxies 
with low SFR or heavily-reddened late-type galaxies with 
a young stellar population and high SFR. 

As developed in Section |5~4} the SFR transition from 
red to blue galaxies can be approximated by a second- 
order polynomial fit in Ms and (U-B) color that is valid 
within the range of SFRs and errors considered in our 
matched DEEP2 / AEGIS sample. The left hand plot of 
Figure [10] shows the calibrated SFRs for the full DEEP2 



sample from 0.74 < z < 1.4 using the [Mb, (U-B), (U- 
B) 2 ] calibration parameter model (see Table p I. Com- 
pared to the S09 SFR contours, the global SFR trends 
are reproduced by this calibration and SFRs smoothly 
transition between galaxies of differing restframe color. 
However, there is some discrepancy for faint blue galax- 
ies; the SFRs become increasingly dependent on (U-B) 
color for (U-B)< 0.5 and M B > -20. This unphysical 
behavior is a consequence of the second-order term in 
the fit where there are few matched S09 data points to 
constrain the SFR calibration. 

The right hand plot of Figure [TU] shows the linear SFR 
calibration using the multicolor calibration with M b , (U- 
B), and (B-V) color used as independent fit parameters. 
Similar to the second-order Mb - (U-B) SFR calibration, 
the multicolor calibration breaks some of the degeneracy 
in intrinsic galaxy color and dust reddening and therefore 
produces a smooth transition in (U-B) color between red 
and blue galaxies. However, the multicolor calibration 
remains insensitive to restframe colors for (U-B)< 1 and 
therefore does not produce an unphysical turnover in the 
SFR dependence for faint blue galaxies. We also note 
that multicolor SFR calibration produces increased SFRs 
for the brightest red galaxies (Mg < —22.5 and (U-B)> 
1). The S09 data shows that bright star-forming galaxies 
with heavily-reddened restframe colors do exist in this 
region. However, because the transition from high to 
low SFR is a strong function of (U-B) color, the SFRs in 
this region will not be highly accurate on an individual 
galaxy basis, but the SFR trend for the entire galaxy 
population will on average be correct. 

While the two calibrations used in Figure [10] are sta- 
tistically equivalent in comparison to the S09 matched 
sample, we opt for the linear multicolor calibration to 
produce be tter qualitative agreement for faint blue galax- 
ies. Figure [TT| shows the SFR produced by the multicolor 
calibration matches the "true" SED SFRs very well, with 
a residual scatter of a a a = 0.28 dex for all galaxies with 
no offset in the mean of the distribution. In fact, the 
scatter is slightly reduced relative to the residual scatter 
seen in the original volume-limited calibration generated 
between 0.74 < z < 1. The additional high rcdshift 
galaxies from 1 < z < 1.4 are predominantly blue galax- 
ies which have a more accurate SFR calibration than red 
galaxies through Mb- The additional blue galaxies re- 
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Fig. 10. — The diagnostic M_e - (U-B) plane with calibrated 0.74 < z < 1.4 DEEP2 SFRs using Mb and a second-order polynomial in 
(U-B) restframe colors (left) and the SFR calibration using Mb, (U-B), and (B-V) restframe colors (right). 
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Fig. 11— SFR fit residuals using MB, (U-B), and (B-V) as 
the independent parameters in the SFR calibration for all DEEP2 
galaxies from 0.74 < z < 1.4. 

duce the average residual scatter for the combined galaxy 
sample. 

4. COMPARISON TO LOCAL [O II] SFR CALIBRATIONS 

The [O II] emission line has long been studi ed as a di- 
agnostic for star formation in active galaxies ( Gallagher 
et~aT[l989; Ho gg et al.|[i"998] [Kelmlc^[T998l |Teplitz et| 
al.| |2003p . To remove the effects of dust extinction and 
oxygen abundance / metallicity, [O II] is typically cali- 
brated to Ha emission in the local universe, where Ha 
luminosity is directly proportional to the number of ion- 
izing photons available from massive stars. The Ha lumi- 
nosity is also more robust against systematic error than 
other emission lines as it is generally insensitive to the 
ISM metallicity, experiences less dust attenuation, and 
absorption due to the underlying stellar continu um can 
be easily modeled ( |Moustakas fc Kennicutt 12006 ). A cal- 
ibration of the [O II] /Ha ratio can be readily performed 
to good accuracy at low redshifts where both lines can 



be simultaneously observed (Kewley et al. 2004| M06). 
At higher redshifts, the Ha line typically becomes diffi- 
cult to measure as it moves to NIR wavelengths, leaving 
estimates of the galaxy SFR dependent on assumptions 
of the mean [O II] /Ha ratio, extinction, and metallicity 
of lower-redshift samples. 

In this section, we study two empirical [O II]-based 
SFR calibrations measured in the local universe and com- 
pare the results t o th e multicolor z ~ 1 SFR calibra- 
tion developed in §3.4| We reiterate that these empirical 
[O II] SFR diagnostics are measured from large galaxy 
samples and are intended to represent the mean trends 
of the galaxy population. Applying the SFR calibrations 
on a individual galaxy-by-galaxy basis may lead to large 
variations and systematic error when compared to indi- 
vidual galaxy SFR diagnostics. Additional measures of 
dust extinction and metallicity should be obtained to ac- 
curately calculate an individual galaxy SFR. 

4.1. [O II] SFR calibration through M B 

To mitigate the systematic effects that limit [O II] 
as a SFR tracer, M06 uses an empirical calibration of 
the [O II] SFR against the restframe .B-band magni- 
tude, which is correlated with both the galaxy extinc- 
tion and metallicity. Mb therefore serves as a proxy for 
the variation in the [O II] /Ha ratio from the local av- 
erage value. The M06 calibration produces good agree- 
ment (0.28 dex RMS scatter) with Ha-based SFRs at 
low redshift. M06 performed a comparison with galax- 
ies at higher redshift (z ~ 0.7 — 1.5) and found that the 
dust-corrected [O II] /Ha ratio remains roughly constant, 
although brighter galaxies tend to have lower ratios (see 
Fig. 22 of M06). The trend of lower observed [O II] /Ha 
ratios for brighter and more massive star-forming galax- 
ies has been further confirmed with NIR measurements 
and tied to increased reddenin g through indepen dent 
measurements of the extinction (Wci ner et al.|[2007 ). To 
obtain SFRs for our sample, we follow the recommen- 
dation of M06 and interpolate between the median SFR 
values (their Table 2) using L[0 II] as a tracer of the SFR 
and Mb to correct the SFR for dust extinction. The ab- 
solute i?-band magnitudes assumed in M06 are in the 
Vega system, but we convert them to AB magnitudes for 
consistency within this study. 

We stress that one should not directly apply the 



11 




-18 



-19 



-20 -21 
M B [mag] 



-22 



-23 



Fig. 12. — The observed [O II] emission line luminosity as 
a function of restframe S-band luminosity for the local sample 
measured in M06 (black triangles) and for DEEP2 galaxies from 
0.74 < z < 1.4 (gray contours). The vertical and horizontal lines 
mark the volume limits for blue DEEP2 galaxies at z = 1 (dot- 
dashed) and z = 1.4 (dotted). The solid (black) line is a weighted 
linear fit to the M06 data in the plotted range, and the dashed line 
is the same fit with a L[0 II] zeropoint shifted by +0.35 dex. 

locally-derived M06 SFR calibration using both observed 
L[0 II] and restframe measured for galaxies at inter- 
mediate redshift. Applying the local SFR relations and 
dust corrections to galaxies at intermediate redshifts can 
introduce considerable systematic uncertainty into the 
SFR calculation. It is well known that the S-band lumi- 



nosity function has evol ved since z = 1 (|Willmer et al. 

f 



2006||Faber et al.||2007h , where the mean M B for inte 
grated galaxy samples has dimmed by Q « 1.3 magni- 
tudes between 2 = 1 and z = 0. To show the discrepancy 
between the M06 local SFR calibration and that of a 
similar L[0 II I an d Mb relation constructed at z ~ 1, we 
plot in FigurejT2"|the raw M06 data for star-forming local 
galaxies (Mb < —18) and 0.74 < z < 1.4 blue galaxies 
from DEEP2. The blue lines in the figure correspond to 
the approximate volume limits of the DEEP2 measure- 
ments at z < 1 (dot-dashed line) and z < 1.4 (dotted 
line). For reference, we plot the weighted linear fit to the 
M06 data set with a black solid line, showing the trend 
of the galaxy continuum luminosity to be correlated with 
the observed [O II] luminosity and therefore SFR for blue 
galaxies. Clearly, the DEEP2 data extends to higher lu- 
minosities in both L[0 II] and Mb than M06, making the 
application of a SFR calibration at intermediate redshift 
heavily dependent on an extrapolation from the local val- 
ues. It is also clear that the DEEP2 blue galaxy samples 
lie well above the mean local trend, however it is not im- 
mediately apparent whether the DEEP 2 luminosities are 
brighter due to an evolution in Mb, L[0 II], or both. 

As sum ing no evolution in M^, the dashed line in 
Fig. [12"1 shows that the aggregate sample of DEEP2 
blue galaxies has a ~ 0.35 dex increase in observed 
[O II] luminosity, and one would incur a similar bias 
in the mean SFR of z ~ 1 galaxies. A more accurate 
calibration of the SFR at z ~ 1 using the M06 data 
should at minimum account for the evolution in galaxy 
luminosity per fixed stellar mass. We therefore correct 
the measured values Mb from DEEP2 galaxies to 



include a dimming factor of Q = 1.3 magnitudes per 
unit redshift before applying the local M06 calibration. 
This simple luminosity correction factor implicitly 
assumes that both the [O II] /Ha ratio and stellar mass 
(or the mass-metallicity relation) has evolved less than 
the galaxy lum inosity, an assumption in agreement with 
recen t studies (Kobulnicky & Kewley 2004 Zahid et al. 
2011 J. Moustakas et. al, 2011, in preparation). 



4.2. [0 II] SFR calibration through Stellar Mass 

Recent efforts to use L[0 II] as a SFR diagnostic have 
focused on incorporating stellar mass into the calibra- 
tion. Stellar mass has been shown to be correlated with 
many galaxy properties, including luminosity, metal- 
licity, dust extinction, and the SFR in SDSS galaxies 



(iBrinchmann et al. 


12004 


Blanton & Moustakas 


2009 


|Garn et al. 


2010). 


The stellar mass estimate serves as 



accounts for the systematic differences in t he observed 



[O II] /Ha ratio. Following the formalism of Gilbank et 



al. (2010 hereafter G10), we co nstruct an empiri cal caT- 



toation of the [O II] SFR from |Kennicutt] ( |1998[ ) using 



SFR 



10 



OAAi 



L[0 II] 



L[0 II] 



[O II} obs /[Ha} obs 1.4 x 10- 41 erg s" 



(3) 



in log(M yr" 1 ), where [O II] f, s /[Ha] & s is the average 
reddened line flux ratio between [O II] and Ha for the 
total sample of galaxies and Aho is the extinction of Ha 
line emission due to dust. In the case of DEEP2 galax- 
ies, we use a ratio of observed [O II]/Ha=0.48 measured 
indirectly through [O II]/H/3 from 0.73 < z < 0.87 and 
Hg/H a from 0.33 < z < 0.39 at M H =-21 (|Weiner et al. 



2007J. The value of [O I I] /Ha is in agreem ent with local 
measurements found in |Kennicutt] ( |1992[ ) , but we cau- 
tion that this value is not singular and can vary greatly 
depending on the galaxy sample selection. For the dust 
attenuation Ah a , we use the mass-de pendent dust ex- 
tinction model of Garn & Best (2010), who performed 
a PCA analysis of SDSS DR7 galaxies using SFR, stel- 
lar mass, and metall icity to produc e a best fit model of 
the Ha absorption. |Garn fc Best| found that a simple 
polynomial model in stellar mass best described the dust 
extinction using 



A 



H, 



0.91 + 0.77X + 0.11A 2 - 0.09A 3 



(4) 



where X = log^/lO 10 M ). Implicitly, this model ac- 
counts for deviations from the observed [O II] /Ha ra- 
tio and the observed differences in metallicity and lu- 
minosity that are correlated with the dust extinction. 
The data used to construct the model cover a large 
dynamic range of stellar mass for "normal" galaxies, 
8.5 < log(M*/M sun ) < 11.5, which spans the majority 
of the z ~ 1 DEEP2 sample. 

4.3. SFR Comparisons 

4.3.1. SED-fit SFRs compared to L[0 II] SFRs 

To understand how the M06 and G10 calibrations com- 
pare to actual SFRs measured at z ~ 1, we again use 
the S09 SFR data from our volume-limited DEEP2 / 
AEGIS sample. Because these calibrations were derived 
from local samples of star-forming galaxies, we restrict 
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Fig. 13. — Comparison of the S09 SFRs to calculated SFRs using 
L[0 II] measurements and empirical dust corrections from the cur- 
rent epoch. The dashed (red) lines show the weighted linear least 
squares fit of the data points and the solid (black) line is the 1:1 
correspondence. The top plot shows the S09 SFRs compared to 
the M06 empirical SFR relation using L[0 II] and M s . The M06 
calibration has been adjusted by Q = 1.3 magnitudes per unit red- 
shift in Mb to account for the luminosity evolution at z ~ 1. The 
dotted (red) line shows the M06 calibration before the luminosity 
evolution is taken into account. The bottom plot shows the S09 
SFRs and the G10 SFRs which uses L[0 II] and stellar mass for 
calibration. 

the sample to only galaxies with blue restframe colors 
and convert all SFRs to a Salpeter IMF for these com- 
parisons. 

Figure [i~3] shows the comparison between the L[0 II]- 
based SFRcalibrations (M06 and G10) and the S09 SFRs 
from 0.74 < z < 1.0. The dashed lines in the figure show 
the linear least squares fit of the sample weighted by the 
estimated errors in both measures of the SFR. The upper 
panel of the plot shows the M06 calibration produces a 
0.29 dex RMS residual scatter relative to the S09 data 
with a mean offset of -0.02 dex. The RMS scatter mea- 
sured from the M06 calibration extrapolated to our sam- 
ple is 0.08 dex larger than the residual error generated 
from the blue galaxy SFR calibration using L[0 II] and 
Mb (see ^3.3). These results indicate that extrapolat- 
ing the SFR behavior through L[0 II] and Mb from the 
local relations to z ~ 1 can contribute to the residual 



scatter, which is not unexpected given the systematic 
uncertainty in the galaxy dust corrections and average 
£?-band luminosity evolution. We also plot a dotted line 
to represent the best fit trend of the M06 calibration 
without accounting for the luminosity evolution between 
z = 0.88 and the current epoch. The difference between 
the SFR distribution means with and without luminosity 
evolution results in a bias of 0.25 dex. 

The bottom panel of Figure 
the G10 calibrated SFRs to t 



13 shows the comparison of 
We find that 



;he S09 SFRs 

the SFR distributions agree with 0.27 dex RMS scatter 
and have a mean offset of -0.06 dex. While the mean 
offset is larger than the M06 calibration, the scatter is 
slightly smaller in the G10 case. Remarkably, the relation 
between stellar mass and Ana in Equation [3] produces a 
relatively accurate calibration at z ~ 1 without having 
to account for any evolution. This result indicates that 
the dust extinction correction required for [O II]-based 
SFRs is better correlated and more constant with stellar 
mass than restframe Mb- 

We note that the slopes of the weighted linear fits 
shown in Figure [l~3| are not well matched to a 1:1 cor- 
respondence. This is due in part to the large amount 
of scatter (~ 0.28 dex) and SFR measurement error rel- 
ative to the overall SFR range probed (~ 1.5 dex) in 
the volume-limited blue galaxy sample. In order to bet- 
ter test the overall accuracy of the [O II] SFR calibration 
slope, we must increase the sample size and probe a wider 
range of SFRs. 

4.3.2. Multicolor SFRs compared to L[0 II] SFRs 
We now compare the multicolor SFR calibration pre- 



sented in £3.4 to the M06 and G10 SFR calibrations ap- 
plied at intermediate redshifts, where we have extended 
the calibrations to include both red and blue galaxies 
from 0.74 < z < 1.4 with measured L[0 II]. As men- 
tioned in |3.1| determining a SFR for a red galaxy based 
on the observed [O II] luminosity can have large sys- 
tematic uncertainty. However, there are some applica- 
tions for which determining a rough SFR for large galaxy 
samples can still produce interesting global measures, e.g 
in measuring how th e SFR correlates with environment 
(Cooper et al. 2006). It is also interesting to check for 
systematic biases m [O II] SFR calibrations when they 
are applied to red galaxies at z ~ 1. 



As shown in the top panel of figure 14 we find that 
the M06 empirical calibration from L[0 II] and Ms is a 
reasonable match to our multicolor SED-fit SFR calibra- 
tion after a c orre ction is applied for luminosity evolution 
in Ms (see £4.1). The residual median offset between 



the SFR calibrations is 0.05 dex with a RMS residual er- 
ror or 0.26 dex after the offset is removed. The residual 
RMS scatter has been reduced slightly over the previ- 
ous volume-limited blue galaxy sample because we have 
added in many star-forming blue galaxies at z > 1. The 
SFR calibration slope differs by 11% per decade, perhaps 
due to additional systematic bias such as an evolution in 
the [O II] /Ha ratio between z — 1 and the current epoch. 
Alternatively, the disagreement could be due to differen- 
tial luminosity evolution as our assumption of a constant 
Q = 1.3 mag/z corresponds primarily to galaxies. 

The bottom panel of Figure [T4| shows that the G10 cal- 
ibration, modified to correct dust extinction in L[0 II] 
through the stellar mass, is an excellent match to our 
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Fig. 14.— Comparison of the M06 and G10 L[0 Il]-based SFR 
calibrations to the multicolor SFR calibration for 0.74 < z < 1.4 
DEEP2 galaxies. The dashed (red) lines show the weighted linear 
regression of the data points and th e solid (black) line is the 1:1 
correspondence. Similar to Figure [T3"| the M06 calibration has been 
adjusted by Q = 1.3 magnitudes per unit redshift in Mj to account 
for the luminosity evolution at z ~ 1. 



multicolor SFR calibration. The mean of the SFR dis- 
tributions have a very small offset of 0.01 dex and the 
slope of the SFRs are in agreement to better than 2% 
per decade. Overall, both the M06 and G10 SFR calibra- 
tions seem to be generally consistent with our multicolor 
SFR calibration, an encouraging result considering the 
multicolor calibration is based only on broadband opti- 
cal colors and restframe magnitude. We do not observe 
a significant systematic bias in the red galaxy SFRs es- 
timated through either L[0 II] -based SFR calibrations 
relative to our multicolor SFR calibration. 

Finally, we compare the SED-fit SFR contours in the 
restframe color-magnitude plane to those generated f rom 
the M06 and G10 calibrations. As shown in Figure 15 
both L[0 II]-based calibrations produce similar behavior 
as the S09 SFR contours in the color-magnitude plane, 
but the G10 calibration is slightly better at reproduc- 
ing the color-independent SFR behavior seen in the blue 
cloud and therefore is in better agreement wit h ou r multi- 
color calibration for blue galaxies (c.f. Figure 14). Com- 



bined with the previous comparison to the SED-fit SFRs, 
these results show that our SFR calibration using only 
restframe optical colors and B-band luminosity are at 
least as accurate as using a local L[0 II]-based SFR re- 
lation and extrapolating it to z ~ 1. 

5. CONCLUSIONS 

We present a calibration of the star formation rate for 
DEEP2 galaxies using SED-fit SFRs from AEGIS and 
observed correlations in [O II], Mb, and restframe col- 
ors. Within a volume-limited sample of galaxies between 
0.74 < z < 1.0, we formulate and minimize a weighted x 2 
fit with respect to these parameters to determine a sim- 
ple, robust empirical SFR calibration. This calibration is 
useful for estimating the global SFR behavior for a large 
sample of optically-measured galaxies at z ~ 1, but it is 
not recommended for use on an individual galaxy basis. 
We summarize our main results as follows. 

1 . Red sequence and blue cloud galaxies at z ~ 1 have 
large differences in their SFRs. As such, an accu- 
rate SFR calibration must minimally separate the 
galaxies into two populations based on restframe 
color or include the restframe color into the cali- 
bration. 

2. We find that a single linear relation in Ms and 
restframe (U-B) and (B-V) colors breaks the ob- 
served degeneracies in the intrinsic galaxy color 
and dust reddening and results in a SFR calibra- 
tion within the accuracy of the available data. The 
calibration produces fit residual errors with 0.27 
dex RMS scatter for the total sample (auue = 0.20 
dex, o~ re d — 0.47 dex) and no significant system- 
atic residual correlations with L[0 II], Mb, or the 
color-M/L stellar mass. Additional constraint from 
observed [O II] luminosity does not improve the 
calibration accuracy, likely due to the fact that the 
SED-fit SFRs are highly correlated with the galaxy 
UV/optical continuum for blue galaxies and L[0 II] 
is less correlated with star formation in red galax- 
ies. 

3. The multicolor SFR calibration developed from the 
matched DEEP2 / AEGIS sample can be extended 
to all DEEP2 galaxies from 0.74 < z < 1.4 without 
a significant increase in SFR residual scatter. The 
fit SFR behavior in the Mb - (U-B) color plane ac- 
curately reflects the average galaxy population in 
the S09 data without introducing significant dis- 
continuities between red and blue galaxies. 

4. We compare the SED-fit SFRs to two calibrations 
(M06 and G10) which use local [O II] luminosity 
measurements as a primary indicator of SFR and 
apply empirical dust corrections through observed 
Mb or stellar mass. We find that the RMS scatter 
between the calculated SFRs are 0.29 dex (M06) 
and 0.27 dex (G10), which is - 0.08 dex larger 
than the residual scatter generated from the best 
SFR calibrations in this study for the same volume- 
limited blue galaxy sample. 

5. The empirical L[0 II] - Mb SFR calibration found 
in M06 should be corrected for the luminosity evo- 
lution in star-forming galaxies when applied at 
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Fig. 15.— The M s - (U-B) plane with calibrated 0.74 < z < 1.4 DEEP2 SFRs using the M06 empirical SFR calibration (left) and the 
G10 empirical SFR calibration (right) based on local SFR measurements. The contour lines correspond to the S09 SFRs of log(i/>)=[0.3, 
0.6, 0.9, 1.2] M yr- 1 . 

z ~ 1. Applying the locally- measured M06 relation United States Department of Energy under contract No. 
at intermediate redshift and not correcting for the DE-AC02-05CH11231. 
mean galaxy luminosity evolution will overestimate 
the z = 1 SFRs by 0.3 dex or more. An L[0 II]- 
based SFR calibration that uses stellar mass to cor- 
rect for the dust extinction (G10) does not require a 
large correction for evolution, indicating that stel- 
lar mass is better correlated with extinction than 
Mb and has less evolution from z ~ 1 to the cur- 
rent epoch. 



6. Finally, we compare the multicolor SFR calibration 
recommended in this study to the L[0 II]-based 
SFR calibrations from M06 and G10. We find 
that our empirical calibration is in good agreement 
with both of the L[0 II] calibrations for all galaxies 
within the SFR calibration accuracy. We conclude 
that the multicolor SFR calibration fit from the 
AEGIS SED-fit SFRs is equivalent in accuracy to 
the local L[0 II] SFR calibrations extrapolated to 
z ~ 1. 
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